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Abstract 



Aims. We present the analysis of 11 nights of V and / time-series observations of the globular cluster NGC 1904 (M 
79). Using this we searched for variable stars in this cluster and attempted to refine the periods of known variables, 
making use of a time baseline spanning almost 8 years. We use our data to derive the metallicity and distance of NGC 
1904. 

Methods. We used difference imaging to reduce our data to obtain high-precision light curves of variable stars. We then 
estimated the cluster parameters by performing a Fourier decomposition of the light curves of RR Lyrae stars for which 
a good period estimate was possible. 

Results. Out of 13 stars previously classified as variables, we confirm that 10 are bona fide variables. We cannot detect 
variability in one other within the precision of our data, while there are two which are saturated in our data frames, 
but we do not find sufficient evidence in the literature to confirm their variability. We also detect a new RR Lyrae 
variable, giving a total number of confirmed variable stars in NGC 1904 of 11. Using the Fourier parameters, we find 
a cluster metallicity [Fe/H]zw = —1.63 ± 0.14, or [Fe/H]uvES = —1.57 ± 0.18, and a distance of 13.3 ± 0.4 kpc (using 
RR0 variables) or 12.9 kpc (using the one RR1 variable in our sample for which Fourier decomposition was possible). 
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1. Introduction 

The ACDM cosmological paradigm predicts that the Milky 
Way formed through t he merger of small gala xies, following 
a hierarchical process ([Diemand et al.l [20071) . As potential 
survivors of these processes, globular clusters are therefore 
important probes into the formation and early evolution 
of the Milky Way. Furthermore, they also allow us to delve 
into the structure of the Galaxy, f or example in con straining 
the shape of the Galactic Halo (|Lux et al.l l2012t) . and to 
study stellar populations. 

There are various methods in use to determine a clus- 
ter's properties, including analysis of its colour-magnitude 
diagram, or obtaining spectroscopy of giant stars in the 
cluster to determine their characteristics. Another way to 
determine the metallicity and distance of a cluster, as well 
as to obtain a lower limit on its age, is to study the pop- 
ulation of RR Lyrae variables that reside within it. By 
obtaining photometry with sufficient time resolution, the 
light curve of the RR Lyrae stars can be analysed with 
Fourier decomposition, yielding parameters that can then 
be used to determine many of the stars' intrinsic proper- 
ties, thanks to various empirical or theoretical relations 
derived from collected observations or theoretical models 



(ISimon fc Clementl [19931 lJurcsik fc Kovacsl 1996; Jurcsik 
119981: iKovacsl 119981: iMorgan et al.ll2007l ). These can also 
then be used as proxies for the properties of the stars' host 
cluster. 

In this paper we use this method to determine the prop- 
erties of the RR Lyrae stars in NGC 1904 (M 79; a = 
05''24 m 10 s ,(5 = -24°31'27" at J2000.0), a globular cluster 
at ~ 13kpc, with metallicity [Fc/H]^ —1.6, and a promi- 
nent blue horizontal branch. We use Fourier decomposi- 
tion to constrain the properties of the cluster itself, follow- 
ing p r evious papers on oth e r clusters (lArellano Ferro et al.l 
20041: lLazaro et al.l 120061; lArellano Ferro et al.l l2008allbl. 



20ld iBramich et aill2011l; lATellano Ferro et alj|20 111) .This 

cluster has been poorly studied with no major published 
long-baseline time-series CCD observations, although it 
contains 13 stars clas sified as variables, in cluding 5 detected 
in a recent study by lAmigo et al.l (|201ll) . Here we present 
an analysis using data spanning almost 8 years, allowing 
us to improve the accuracy of the RR Lyrae periods and 
to detect variations in period or amplitude that might be 
indicative of the presence of the Blazhko effect. 

Although it has been suggested N GC 1904 belonged t o 



nkains@eso.org 



the Canis Major ( CMa) dwarf galaxy (jMartin et al.l l 
other studies (e.g. iMateu et al.l 120091 ) have questioned this 
due to the lack of blue horizontal branch stars in the CMa 



1 



N. Kains et al.: Constraining the parameters of globular cluster NGC 1904 from its variable star population 



Date 


AT 

ivy 


ty (sj 


AT 
JVJ 


ti (s) 


l W W L 1 1 f'\f~\ 1 

20041004 


7 


60-200 







20041005 


3 


100 







20090108 


16 


150 


14 


120 


20100123 


19 


80-150 


21 


40-120 


20100221 


16 


100-200 


14 


50-100 


zUlUUoU / 


10 


on 1 on 
yU-lzU 




40-00 


20110312 


16 


ion 1 /i n 
120-140 







20110313 


10 


80-150 







20111104 


30 


90-300 


32 


25-120 


20120301 


10 


150-180 


10 


50-65 


20120302 


5 


100-120 


5 


45 


Total: 


147 




109 





Table 1. The distribution of observations of NGC 1904 
for each filter, where the columns Ny and Nr represent the 
number of images taken for the filters V and /, respectively. 
We also provide the exposure time, or range of exposure 
times, employed during each night for each filter in the 
columns ty and t\. 



colour-magnitude diagram (CMD) compared to this and 
other clusters also potentially associated with the CMa 
overdensity, such as NGC 1851, NGC 2298 and NGC 2808. 
Further study of these clusters is therefore crucial to im- 
proving our understanding o f the structure of the Galactic 
disc, as iMartin et al.l (|2004h also suggested that CMa is 
made up of a mixture of thin and thick disc and spiral arm 
populations of the Milky Way. 

In Sec. [2J we detail our observations, and in Sec. [31 we 
discuss variable objects in NGC 1904 both in the literature 
and those which we detect in our data. In Sec. [4] we perform 
Fourier decomposition of the light curves of some of the 
variables to determine their properties. We then use these 
in Sec. [5] to calculate values for the cluster parameters. We 
place these values into context in Sec. H] by comparing them 
to those found for other clusters in the literature. Finally, 
we also discuss briefly the peculiar spatial distribution of 
the RR Lyrae variables that we find in NGC 1904. 

2. Observations and reductions 

2.1. Observations 

We obtained Johnson F-band data of NGC 1904 using the 
2.0m telescope at the Indian Astrophysical Observatory in 
Hanle, India, on 11 nights spanning from April 2004 to 
March 2012, a baseline of almost 8 years. The cluster was 
also observed in the / band, except in 2004 when no / data 
was taken. This resulted in 147 V^-band and 109 /-band 
images. The observations are summarised in Table [TJ 

The CCD used was a Thompson 2048 x 2048 pixel, with 
a field of view of 10.1 x 10.1 arcmin 2 , giving a pixel scale 
of 0.296 arcsec per pixel. Given the distance to the cluster 
we derive later in this paper, this corresponds rougly to an 
area of 40 pc x 40 pc centred on the core of the cluster. 

2.2. Difference Image Analysis 

We reduced our obs ervations using the DanDl43 pipeline 
(|Bramichll2008l l2012t ) to obtain high-precision photometry 

1 DanDIA is built from the DanlDL library of IDL routines avail- 
able at |http : / / www . danidl . co . uk | 
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Figure 1. Linear regression fits to the relation between in- 
strumental and Johnson-Kron-Cousins magnitudes for the 
V (top) and / (bottom) bands. 



of sources in the images of NGC 1 904, as was done in pre- 
vious globular cluster studies (e.g. iBramich et al.|[2~012l see 
that paper for a detailed description of the software). We 
produce a stacked reference image in each band by select- 
ing the best-seeing images (within 10% of the best-seeing 
value) and taking care to minimise the number of saturated 
stars. Our resulting reference image in the V filter consists 
of a single image with an exposure time of 100 seconds, 
and full width half-maximum (FWHM) of the point spread 
function (PSF) of 5.37 pixels (~ 1.5"). In the / band, the 
reference image is made from 3 stacked images with a total 
exposure time of 150 seconds and a PSF FWHM of 4.79 
pixels (~ 1.4"). 

Because of the way difference imaging works, the mea- 
sured reference flux of a star on the reference image 
might be systematically too large due to contamination 
from other nearby objects (blending). Since non- variable 
sources are fully subtracted on the difference images, how- 
ever, this problem does not occur with the difference im- 
ages. Hence variable sources with overestimated reference 
fluxes will have an underestimated variation amplitude, al- 
though crucially the shape of the light curve is unaffected. 
Furthermore, some of the targets which are in the vicin- 
ity of saturated stars or are saturated themselves will have 
less precise or non-existent photometry. Although we tried 
to minimise the number of saturated stars in our reference 
images, some remained and this affected the photometry of 
some of the objects in this cluster. 
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# Filter HJD M std m ins a m /ref cfref /diff Cdiff P 

(d) (mag) (mag) (mag) (ADU s" 1 ) (ADU s" 1 ) (ADU s' 1 ) (ADU s" 1 ) 

V3 V 2453283.39244 16.235 16.875 0.005 2531.207 6.272 -740.195 7.926 0.9831 

V3 V 2453283.39990 16.237 16.877 0.003 2531.207 6.272 -748.391 4.372 0.9900 

V3 I 2454840.11764 15.555 17.093 0.005 1206.297 6.739 380.674 9.788 1.5286 

V3 I 2454840.12374 15.566 17.104 0.006 1206.297 6.739 358.597 11.513 1.5294 



Table 2. Format for our time-series photometry, for all confirmed variables in our field of view. The standard M st d and 
instrumental mi ns magnitudes listed in column 4 and 5 respectively correspond to the variable star, filter and epoch of 
mid-exposure listed in columns 1-3, respectively. The uncertainty on mj ns is listed in column 6, which also corresponds to 
the uncertainty on M st d- For completeness, we also list the reference flux / re f and the differential flux /diff (columns 7 and 
9 respectively), along with their uncertai nties (columns 8 and 10), as well as the photometric scale factor p. Definitions 
for these quantities can be found in e.g. iBramich et al.l (|201ll ). Eq. 2-3. This is a representative extract from the full 
table, which is available with the electronic version of the article (see Supporting Information). 



2.3. Photometric Calibrations 

Instrumental magnitudes were converted to standard 
Johnson-Kron-Cousins magnitudes by fitting a linear re- 
lation betw een the known magnitudes of standard stars in 
NGC 1904 (|Stetsodl2000l ) and the light curve mean mag- 
nitudes. These relations, plotted in Fig. [TJ were then used 
to convert all instrumental V— and I— band light curves 
to standard magnitudes. We note that the standard stars 
we use cover the full magnitude and colour ranges of our 
CMD, justifying our choice not to include a colour term in 
our calibration relations. Calibrated data for all variables 
is available with the electronic version of this paper, in the 
format given in Table O 



2.4. Astrometry 

We used the online tool astrometry.net (| Land [200 9t ) to 
obtain an astrometric fit to our V— band reference im- 
age, and used this to calculate the J2000.0 coordinates of 
all objects we discuss in this paper. This tool uses the 
USNO-B catalog of astrometric standards to find an as- 
trometric fit to images. The coordinates are those of the 
epoch of our V— band reference image, which was taken at 
HJD=2453284.41 d. We estimate that the astrometric fit is 
accurate to within ~ 0.1" in a and ~ 0.3" in S. 



3. Variables in NGC 1904 

Th e first five var iables (Vl-5) in this cluster were reported 
bv iBailevi (Il902l) . examining plates taken at the Harvard 
College station in Arequipa, Peru. Of these, VI has not 
been confirmed as variable by subseque nt studies, wh ile V2 
is classified as a semi-regular variable bv lRosin ol ffl952h. One 
more variable (V6) was then identified by Rosinol ( |1952l ) 



from observations carried out between 1948 and 1951 with 
the 24-inch reflector of Lojano; these also provided period 
estimates for the new variable as well as two of Bailey's 
original variables (V3-4). Rosino also corrected the coordi- 
nates for V3 and V4, and commented that V5 in Bailey's 
original study was not detectable in the new observations 
due to how crowded his photographic plates were. 

Two mo r e var iables (V7-8) were reported by 
Sawyer Hoed (119731) . although the online catalogue of 
Clement et al.l ( 200lh notes that these are unlikely to be 



RR Lyrae variables, as they sh ould have been dete cted as 
such in the subsequent work of lAmigo et al.l (|2011l ). That 
study analysed -B-filter observations taken in 2001 at the 
Danish 1.54m telescope in La Silla using image subtraction, 
and identified an additional five new RR Lyrae variables 
in NGC 1904, V9-V13. They also estimated the periods 
for these, as well as the periods of V3, V4, V5 (which they 
mistakenly referred to as "NV6") and V6. 

This amounts to a total of 13 stars classified as vari- 
ables in this cluster, 9 of which have had their period es- 
timated, and 3 of which are suspected variables. Only the 
most rec ent time-series stu dy was carried out using CCD 
cameras (jAmigo et al.ll20fl"l ) , and no work employing multi- 
band time series observations for this cluster with CCD 
cameras has been published so far. 

In this paper we use the notation introduced by 
lAlcock et al.l (|2000f ) , referring to fundamental-mode pulsa- 
tion RR Lyrae stars as RR0 and to first-overtone RR Lyrae 
pulsators as RR1, rather than RRa6 and RRc, respectively. 



3.1. Stars that do not show variablity 

VI did not show signs of variability in our data, to within 
the limit set by its V^-band root mean square (rms) scat- 
ter of 0.015 mag and /-band r ms of . 044 m ag. This is in 
agreement with the finding of iRosin" 3 (119521V which men- 
tions that VI might not be variable. 



3.2. Variable candidates without light curves 

V2 and V8 are both saturated in our reference images, and 
we are therefore unable to present light curves for them. To 
the best of our knowledge, there is no published light curve 
for either of these two objects. V8 was only reported as 
variable in a private communication by Tsoo Yu-Hua to H . 



Sawyer Hogg (see online catalogue of 



Clement et al.ll200lD . 
and was not detected by lAmigo et al. (|2011 ) , which means 
it i s unlike l y to b e an RR Lyrae variable. V2 is reported 
by iRosinol (|1952f) as a semi-regular or irregular variable, 
with no period estimate. We suggest that more time-series 
observations are needed to confirm the true variability of 
V2 and V8 and we believe that they should remain variable 
candidates. 
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V9 V10 V11 V12 V13 




Figure 2. Finding chart for the confirmed variable objects in NGC 1904, using our V reference image. North is up and 
East is to the right. The image size is 4.44 x 5.18 arcmin 2 , while each stamp size is 23.7 x 23.7 arcsec 2 . White circles and 
labels indicate the locations of the variables, and each of the variables we detect in our data is shown in an individual 
stamp, its location marked by a cross-hair. V2 and V8 are saturated in our reference image and therefo re we do not 
presen t light curves for them, but their locations are sh own on this finding chart using the coordinates of ISamus et al.l 
(|2009ft . The tidal radius r t ~ 500" (|Lanzoni et al.ll2007D lies outside of this plot, 
4 
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Table 3. Data for all confirmed variable stars in NGC 1904. Periods are determined using phase dispersion minimisation 
(PDM) and the string length method (see text). < V > and < I > are intensity- weighted mean magnitudes. a poten- 
tially RRd, see text. 6 error-weighted mean magnitudes, rather than intensity-weighted magnitudes. Coordinates for all 
confirmed variables are given in Table 01 
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Table 4. J2000.0 equatorial coordinates for all confirmed 
variable stars in NGC 1904, derived from our astromctric 
fit to the V — band reference image. 



3.3. Detection of known variables 

Tablc[3]lists all confirmed variables in NGC 1904, including 
the 10 already known variables, and their coordinates ob- 
tained from the astrometric fit to our reference image are 
given in Table A finding chart of the cluster showing the 
location of the variables is shown in Fig. [21 When possi- 
ble, the variable type is also given in the last column of 
the table. We performed a peri od search using bo th phase 
dispersion minimisation (PDM, Stellingwerf 1978f) and the 
"string length" method (jLafler fc Kinmanlll965t ). The pe- 
riods we find offer an impr ovement in p r ecisio n compared 
to the periods reported bv lAmigo et al.l (|2011), thanks to 



our lo ng baseline. Indeed, using the periods of lAmigo et al.l 
(pOll with our data results in poorly phased light curves. 

V7 shows some long-term variability and no short-term 
variability, and is therefore unlikely to be an RR Lyrae 
variable ; furthermore , it w as not detected as such by the 
work of lAmigo et al.l (|201l[ ). and its position on our CMD 
(Fig. [3]) also makes it an unlikely RR Lyrae candidate. A 
plot of average nightly magnitude against time for V7 is 
shown in Fig. 01 showing the long-term variation. 

Fig. [5] shows phased V light curves of all confirmed vari- 
ables for which a period could be estimated, /-band light 
curves were also obtained, except for V12 as the photom- 
etry was not good enough, with poor seeing and the star 
being located in the centre of the cluster, and V13 whose 
location next to a bright star and poor seeing prevented us 



3500 4000 4500 5000 
HJD-2450000 



3500 6000 



Figure 4. Average F-band magnitude of V7 for each ob- 
serving night, showing some long-term variability of an un- 
determined nature. The X^-band data is plotted in black, 
while the /-band data is plotted in red, with a fixed offset 
to bring the two data sets onto the same plot. The varia- 
tions in (V — I) therefore appear to indicate a change with 
time in the colour of V7, although the /—band data is poor. 
The x-axis is given in HJD-2450000. 



from extracting a good quality /-band light curve. We do 
not show our /—band light curves here, because they are 
noisy and were not significant to the analysis presented in 
this paper. However, the light curves in both V and / are 
available in electronic format. 

The phased light curves of both V 4 and V6 sug gest that 
they might exhibit the Blazhko effect (|Blazkoll907f) , leading 
to variations in the period and amplitude of the variation in 
brightness. The effect is weak in V4, but very pronounced 
in t he case of V6 , and we suggest that this is the reason 
why lAmigo et al.l (|2011l) also had issues phasing the light 
curve for this object. Unfortunately our time coverage docs 
not allow us to estimate the period and magnitude of the 
Blazhko effect for any objects in this study, but we note 
the similarity of the light curves of V4 and V6 to those of 
many of the RR1 variables with a detect ed Blazhko effect 
discovered bv lArellano Ferro et al.l (|2012l ). 

We see from the reference image residuals after PSF 
fitting that V5 and V10 are both blended in our data with 
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Figure 3. V — I, V colour-magnitude diagram of NGC1904 extracted from the reference images. The locations of the 
confirmed RR Lyrae variables with V and I photometry are marked with blue filled circles (RRO variables) and blue filled 
triangles (RR1 variables), V7 with a green filled circle, and VI, for which we do not find evidence of variability, with a red 
cross. Note that V5 and V10 lie above the horizontal branch because they are significantly blended, as discussed in the 
text. Error bars for various mag nitude levels are shown near the right edge of the plot. Also shown are an interpolation 
of the IVandenBerg et al.l (|2006l ) isochrones for a metallicity of [Fe/H]zw = —1-63 (see Sec. 15. 2p at ages of 9, 10, 11, 12, 
13, 14 and 15 Gyr (marked by lines of different colours; see Sec. 15 3[) . 



stars of similar brightness, mostly due to the poor seeing of 
~ 1.5" in our data. This means that although the shape of 
the light curves is not affected (see Sec. I2.2p . the reference 
flux is overestimated for these objects; their position on the 
CMD is also affected, as can be seen in Fig. ■ This is taken 
into account in the analysis when deriving properties of the 
variable stars. 

However, we note here that the position of V5 and V10 
on our CMD could also be due to NGC 1904 having a 
prominent blue HB morphology. Due to this, some of the 
RR Lyrae in this cluster could be the progeny of blue HB 
stars during their evolution to the asymptotic giant branch 
(AGB) phase. Such stars cross the instability strip at higher 
luminosities and could therefore appear on our CMD where 
we find V5 and V1 0. The CMD obtained for this cluster by 
iPiotto et al.1 (|2002t ) indeed shows such stars present in NGC 
1904. However, we also examined HST archival data on this 
cluster and verified that these two stars are indeed blended 
within the FWHM of our reference images. 

Finally, V9 is also not well phased with our best period 
esti mate, and from the shape of its light curve we suggest, 
like lAmigo et al.l ([201 1[ ) , that it could be a double-period 
(RRd) variable. 



3.4. Detection of new variables 

In order to detect new variables, we employed three meth- 
ods: firstly, we constructed a stacked image S consisting of 
the sum of the absolute values of the deviations D of each 
image from the convolved reference image, divided by the 
pixel uncertainty er, so that 



Si 



E 



\D 



kij | 



(1) 



Stars which consistently deviate from the reference im- 
age then stand out on this stacked image. Using this 
method, we discovered V14, a previously unpublished RR1 
variable (see Table [3]). Secondly, we inspected the light 
curves of objects which stood out on a plot of root mean 
square magnitude deviation versus mean magnitude, shown 
in Fig. [51 

Finally, we also searched fo r variables by com puting the 
string-length Sq statistic (e.g. lDworetskvlll983l) for all our 
light curves. We inspected visually all the light curves of 
stars with Sq < 0.3, where the threshold value of 0.3 was 
chosen by inspecting the distribution of Sq (see Fig. [7]). 
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Figure 5. Phased V-band light curves for the variables with a period estimate. Different colours are used for different 
nights, with the colour coding given in the form of a colour bar spanning the whole time range of the observations (top 
panel). Also shown is the Fourier decomposition fit to the lightcurves, for those stars for which this fit was performed 
(see text). The median data error bar is plotted in the top left corner of each plot. 



However, this did not reveal any additional RR Lyrae or 
other variable candidates. 

We estimated the period of V14 using both PDM and 
the string-length method, but like the previously known 
variables V4 and V6, we suggest that V14 is affected by 



the Blazhko effect, which explains the difficulty in phasing 
our light curve. 
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Figure 6. Plot of root mean square magnitude deviation 
versus mean magnitude for all stars for which photometry 
was obtained. RR Lyrae variables are plotted as blue filled 
circles, V7 as a green filled circle, and VI, for which we 
do not find evidence of variability, is shown as a red cross. 
Plots are for the V— band (top) and /—band (bottom). 
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Figure 7. Distribution of the Sq statistic (see text) applied 
to the y-band light curves of all stars in our sample (black 
filled circles). A blue horizontal line marks the threshold 
of Sq = 0.3 below which we inspected the light curves for 
periodic variability. The locations of the confirmed variables 
in our sample are marked with red filled circles; V7 is not 
shown as it does not appear to be periodic, but rather shows 
long-term variability. 



where m(t) is the magnitude at time t, N is the number of 
harmonics used in the fit, P is the period of the variable, E 
is the epoch, and Ak and <pk are the amplitude and phase of 
the k th harmonic. Epoch-independent Fourier parameters 
arc then defined as 



4. Variable Properties 

A (V — I, V) colour-magnitude diagram of NGC1904 ex- 
tracted from the reference images is shown on Fig. [31 with 
the location of the variables marked. Although in theory a 
deeper CMD could be obtained by combining many images, 
in practice the large variations in seeing and the number 
of images with significantly worse seeing than the reference 
image makes this impractical in this case. Variable are plot- 
ted on Fig. [3] using their error- weighted mean magnitudes. 



4.1. Fourier decomposition 

For variables with good enough phase coverage, using 
Fourier decomposition of their light curves can help us de- 
rive several properties of both the variables themselves and 
of t he cluster in which they are located. We exclude V6 from 
the analysis that follows because of the large dispersion in 
our best phased light curve (possibly due to the Blazhko 
effect), but retain V4 and V14, for which the effect is weak. 
Fourier decomposition amounts to fitting light curves with 
the Fourier series 



JY 



i(t) = A + A k cos 



fc=i 



2irk 



'-(t-E) + cf> k 



(2) 



Rij — Ai I Aj 

<t>ij = j4>i - i( t>j 



(3) 
(4) 



We chose to fit the minimum number of harmonics that 
yielded a good fit, taking care not to over-fit light curve 
features. As a check of the dependence on N of the pa- 
rameters we derive for each variable, we also derived stellar 
parameters using Fourier parameters obtained with more 
harmonics, and compared them to those given in the text. 
We found very little variation, with any changes smaller 
than the error bars quoted. 

We list the coefficients Ak we obtain for the first four 
harmonics in Tabic [5J as well as the Fourier parameters 
021 j <^3i and 04i for the variables for which we could ob- 
tain a Fourier decomposit ion. We also list t he de viation 
parameter D m , defined bv I Jurcsik fc Kovacsl (|l996t) to as- 
sess whether fit parameters can be used to derive proper- 
ties of the RR Lyrae variables. Although I Jurcsik fc Kov&csl 
(119961 ) used a criterion whereby fits should have D m < 3 
for their empirical relations to yield reliable estimates of 
stellar properties, a less stringe nt criterion of D m < 5 has 
been used by other authors (e.g. lCacciari et aT1l2005l ). Here 
we also adopt D m < 5 as a selection criterion to estimate 
stellar properties, which excludes one variable, V4, from 
our calculations of the cluster's parameters (Sec. [5]). 



N. Kains et al.: Constraining the parameters of globular cluster NGC 1904 from its variable star population 



# 


A 


A t 


A 2 


^3 


A 4 


021 


031 


041 


N 


D m 


RRO 


V3 


16.049(1) 


0.330(1) 


0.165(1) 


0.107(1) 


0.057(1) 


4.222(9) 


8.680(14) 


6.805(19) 


6 


2.2 


V4 


16.111(1) 


0.408(1) 


0.195(1) 


0.134(1) 


0.062(1) 


4.031(5) 


8.318(8) 


6.599(13) 


7 


7.0 


V5 


15.724(1) 


0.224(2) 


0.114(2) 


0.085(1) 


0.052(1) 


4.115(23) 


8.535(28) 


6.466(32) 


6 


3.5 


V10 


15.027(1) 


0.102(1) 


0.073(1) 


0.018(1) 


0.019(1) 


4.045(16) 


8.644(35) 


6.543(37) 


5 


3.7 


V13 


16.061(1) 


0.273(1) 


0.132(1) 


0.075(1) 


0.020(1) 


4.317(15) 


8.825(22) 


7.290(49) 


4 


2.3 


RR1 


V14 


16.154(1) 


0.137(1) 


0.007(1) 


0.007(1) 


0.005(1) 


4.841(140) 


3.079(121) 


2.226(152) 


4 





Table 5. Coefficients Ak and selected Fourier parameters for the Fourier decomposition fit of the RR Lyrae variables in 
NGC1904. The number of fitted harmonics N is also given. Numbers in parentheses are the 1-a uncertainties on the last 
decimal place. Note that Aq for V5 and V10 is too bright and Ai,A2, A3 and A4 are underestimated because these two 
objects are significantly blended, so we exclude them from calculations in w hich this parameter is us ed to estimate the 
distance to the cluster. We also list the deviation parameter D m defined bv I Jurcsik fc Kovac s (1996) to assess whether 
fit parameters can be used to derive properties of the RR Lyrae variables. 



4.2. Metallicity 
We use the empirical relations of I Jurcsik" 



Kovacsl ((19961) 



to derive the metallicity [Fc/H] for each of the variables for 
which we could obtain a successful Fourier decomposition. 
The relation is derived from the spectroscopic metallicity 
measurement of field RRO variables, and it relates [Fc/H] 
to the period P and the Fourier parameter 0^, where s 
denotes a parameter obtained by fitting a sine series rather 
than the cosine series we fit with Eq. ((2). [Fe/H] is then 
expressed as 



[Fc/H]j = -5.038- 5.394 P + 1.345, 



'31 



(5) 



where the subscript J denotes a non-calibrated metallicity, 
the period P is in days, and 0f ■ can be calculated via 



# 


[Fc/Hjzw 


M v 


log(L/L ) 


T e g 


RRO 


V3 


-1.70(2) 


0.366(3) 


1.288(2) 


6232(7) 


V4 


-1.66(2) 


0.420(2) 


1.291(1) 


6354(4) 


V5 


-1.59(4) 


0.548(4) 


1.228(2) 


6292(12) 


V10 


-1.71(4) 


0.564(2) 


1.195(1) 


6161(14) 


V13 


-1.39(3) 


0.458(2) 


1.262(2) 


6287(11) 


RR1 


V14 


-1.73(6) 


0.584(8) 


1.409(4) 


7293(12) 



Table 6. Physical parameters for the RR Lyrae variables 
calculated using the Fourier decomposition parameters and 
the relations given in the text. Numbers in parentheses are 
the l-o" uncertainties on the last decimal place. Note that 
My for V5 and V10 is overestimated because these two 
objects are significantly blended, so we exclude them from 
the calculation for the distance to the cluster. 



J) 



(6) 



We tran s form these to the metallicity scale of 
Zinn fc Westl (fl98l (hereafter ZW) using the relation from 
Jurcsik! (119951 ): 



[Fe/H] 



zw 



[Fc/H] j - Qj 
1.431 



(7) 



For our only RR1 variable with a good Fourier decom- 
position, V14, we calculated the metal licity using the em- 
pirical relation of iMorgan et al.l <|2007l ). linking [Fe/H], P 
and 



031: 



[Fe/H] zw = 2.424- 30.075 P + 52.466 P 2 (8) 
+0.982 03i + 0.1310^ - 4.198 3 iP 

Metallicity values calculated using Eq. ([5]) & (J7J) and (JSj) 
are given in Table [5] 



4.3. Effective Temperature 

Another property of the RR Lyrae that can be estimated 
thro ugh the Fourie r parameters is the effective tempera- 
ture. [Jurcsik (199<§) derived empirical relations linking the 
(V — K)o colour to P as well as several of the Fourier co- 
efficients and parameters: 



(V-K)q = 1.585 + 1.257 P- 0.273 Ai -0.2340^ (9) 
+0.062 0^ 

logoff = 3.9291 - 0.1112 (V - K) (10) 
-0.0032 [Fe/H] j . 

For RR1 variables, we use a corresponding relat i on de- 
rived from theoretical models by ISimon fc Clement] (|1993|) 
to calculate T g, 



log T cff = 3.7746 - 0.1452 log P + 0.0056 < 



•31 ■ 



(11) 



The temperatures we derive using these relations are 
given in Table [51 note that because V5 and V10 are signifi- 
cantly blended, the temperatures we calculated for them 
are underestimated. As noted by Bramic h et al.l ([201 ll ). 
there are several caveats to deriving temperatures with 
Eq. mnjl and (HU. The values of T oS derived for RRO and 
RR1 variables are on different absolute scales. Furthermore, 
the effective temperatures we derive here deviate sys- 
tematically fi^mthe relations predicted by evolutionary 
models of Casteih] (119991) o r the temperature scales of 
iSekieuchi fc Fukueital (l2000h: however this was already 
noted in previous w ork ( Arellano Ferro et al.ll2008al l2010t 
iBramich et al.ll2~01l[ ). and the temperatures given here are a 
useful comparison with previous analogous studies on other 
clusters. 
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4.4. Absolute Magnitude 

We use the empirical relations of iKovacs fc Walked (pool 
to derive V-band absolute magnitudes for the RRO vari- 
ables. The relation links the magnitude to Fourier coeffi- 
cients through 



-1.876 logP - 1.158 A x + 0.821 A 3 + K , (12) 



where Kq is a constant. A s in iBramich et al.1 (|2011l ) and 
lArellano Ferro et~a!1 (|2010l) , we choose a value of K = 0.41 
mag to be consi stent with a true LMC distance modulus of 
/io = 18.5 mag (iFreedman et al .1 20011) . For RR1 variables, 
we use the relation of Kovacsl ( 1998h . 



My 



-0.961 P- 0.044, 



4.447^4 + -^! 



(13) 



where K\ is a constant, for which we choose a value of 1.061 
with the same justification as for our choice of Kq. 

Wc also converted the magnitudes to luminosities using 

log (L/Lq) = -0.4 [M v + B c (T cS ) - M bol)0 ] , (14) 

where we use M\> \ q = 4.75, and Pc(Pcff) is a bolomctric 
correctio n which we determ i ne by interpolating from the 
values of iMontegriffo et ah! (|1998l ) and using our best-fit 
value of T c ff. We report values of My and log(L / Lq) for 
the RRO and RR1 variables in Table [SJ and also note that, 
because V5 and V10 are heavily blended, their derived value 
of My is likely to be overestimated; we therefore do not use 
them to estimate the distance to the cluster. 



5. Cluster Properties 

5.1. Oosterhoff Type 

Although the number of variables for which the period 
could be determined is small, we calculate the mean periods 
of the RRO and RR1 variables as < P RR0 >= 0.71 ± 0.07 
and < P Rm >= 0.33±0.02. RR1 variables account for 40% 
of all RR Lyrae stars in this cluster. 

The value of < Prro > and the high fraction of RR1 
variables makes this cluster an Oosterh off type II, in agree - 
ment with previous classifications (e.g. lAmigo et al]|2011l) . 
We also show the distribution of variables in a Bailey di- 
agram on Fig. [8l and confirm that the distribution is con- 
sistent with an Oos terhoff type II cluster; as found by 
ICacciari et alJ (|2005l ). type II clusters have log P — Ay 
tracks that correspond to the tracks of evolved stars in 
the "prototype" type I cluster, M3. Although the metal- 
licity value we derive for NGC 1904 makes it o ne of the 
most metal-rich Oosterhoff type II clusters (e.g. ISandagel 
|1993| ). there are several other known type II Galactic clus- 
ters with simila r metallicities such as M2 (|Lee fc Carnevl 
19991) and M9 (IClement fc Sheltonl Il999al). Further more, 



Lee fc Carnevl (|1999f ) and lClement fc Sheltonl (|1999bf) both 



concluded that the Oosterhoff dichotomy was due to evo- 
lution rather than metal content; we therefore do not con- 
clude that our metallicity value and Oosterhoff classifica- 
tio n are contradictory. 

Ivan Albada fc Bakerl (|l973h suggested that in 
Oosterhoff type II clusters, RRc evolve to lower tem- 
peratures, but that there is a hysteresis effect so that 



i .2 ; 
1.0 - 



< 0.6 





0.0 

-0.6 



-0.5 -0.4 -0.3 -0.2 -0.1 -0.0 

log P 



Figure 8. Bailey diagram of the V-band light curve ampli- 
tude versus the logarithm of the period for RR Lyrae stars 
in NGC 1904. RRO variables are plotted as filled circles, and 
RR1 variables as open triangles; V9 is plotted as an open 
square as it is a potential double-mode RR Lyrae star. We 
also mark the locations of V5 and V10, whose amplitude is 
underestimated because of blending. Al s o plo tted as solid 
lines are the relations of ICacciari et al.1 (|2005l ) for evolved 
stars; these are obtained by applying a shift of +0.06 to 
log P in the relations derived for Oosterhoff type I cluster 
M3. 



mode switching occurs at lower temperatures for stars 
evolving to lower temperatures than for stars evolving 
to higher temperatures. Indeed for NGC 1904, we find 
a mean temperature for RRO stars (excluding V4, V5 
and V10) of 6260 ± 39K, which is lower than the means 
found for O osterhoff type I clus ters, e.g. 6418 ± 10K for 



NGC 6981 (IBramich et al 



1997) 



6633 ± 257K for NGC 
6465 ± 22K for M 5 



4147 (lArellano Ferro et al- L 
dKaluznv et al J 1200(1) . 6494 ± 35K for NGC 1851 (IWalker 



I199S ). 6619 ± 64K for NGC 6171 (IClement fc Shelton 



5.2. Metallicity 



Although ISandagel (|2006l ) found that there is a correlation 
between the mean period of the RRO variables and the clus- 
ter metallicity for Oosterhoff type I clusters, this relation 
breaks down when < Prro > is larger than ~ 0.6, i.e. for 
Oosterhoff type II clusters ([Clement et al.ll200ll ). 

On the other hand, we can use the metallicities derived 
in Sec. 14.21 to estimate the cluster's metallicity. Taking an 
average of the RR Lyrae metallic i ties l isted in Table ® 
and assuming as in IBramich et al.1 (|201ll ) that there is no 
systematic offset between metallicity estimates for the dif- 
ferent types of variables, we find a mean metallicity of 
[Fe/H]zw = — 1.63 ± 0.l43- This is in good agreement with 
values in the literature, listed in Table [7] Note that V4 was 
excluded from this calculation as it has D m > 5, meaning 
that its metallicity estimate might be unreliable. The error 
bars given here and for all cluster properties are the scatter 
around the mean value. 



2 Note that the values in Table [6] are rounded off, which ex- 
plains why the mean of those values is not equal to the value we 
derive for the cluster metallicity. 
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We also transform this metallicity value to the scale of 
ICarretta et al.l (|2009(l who derived a new metallicity scale 
based on GIRAFFE and UVES spectra of red giant branch 
(RGB) stars in 19 globular clusters. The transfo rmation 
from the ZW to the UVES (|Carretta et al.ll2009t) scale is 
given as 



[Fe/H] UVES — -0.413 + 0.130 [Fe/H] zw - 0.356 [Fe/H] zw . 

(15) 

Using this we find a metallicity for NGC 1904 of 
[Fe/H]uvES = —1-57 ± 0.18, in e xcellent agreement w ith 
the value found for this cluster by ICarretta et al.l (|2009() of 
[Fe/H] UVES = -1.58 ±0.03. 



5.3. Distance 

Having derived the absolute magnitudes of the RR Lyrae 
stars in Sec. 14.41 we can use the Aq parameter, which cor- 
responds to their mean apparent V magnitude, to derive a 
distance modulus to NGC 1904. Note that My for V5 and 
V10 is overestimated because these two objects are signif- 
icantly blended, so we exclude them from this calculation, 
as we do with V4 because of its value of D m above our se- 
lection threshold. The mean value of the apparent V band 
magnitudes for the RR0 variables is 16.055 ± 0.008 mag, 
while the mean of the absolute magnitudes is 0.412 ± 0.065 
mag. This yields a distance modulus of /i = 15.643 ± 0.066. 
Using the parameters for our lone RR1 variable (V14, see 
Ta bles EfcEt, we find u = 15.570. 

iStetson fc Harris! (|1977t ) used various methods to esti- 
mate the reddening towards this clus ter and obtained val - 
ues of E(B - V) lower than 0.03 mag. iFerraro et all (|l992t) 
found a red dening for NGC 1904 o f E(B - V) = 0.01 ± 0.01 
mag, while iKravtsov et al.l ()1997l ) estimate the reddening 
for this cluster to be E(B - V) = 0.035 to 0.05 mag, de- 
pending on the metallicity of their best-fit isochrone. 

Adopting a value of E(B - V) = 0.01 ± 0.01 mag and 
a value for Ry =3.1 for our Galaxy, we find a mean true 
distance moduli of /in = 15.612±0.066 mag and /io = 15.539 
mag for RR0 and RR1 variables in our samples. This yields 
mean distances of 13.26 ± 0.41 kpc and 12.87 kpc. These 
values are sensitive to the adopted value of E{B — V), and 
using the largest value found in the literature, E(B — V) = 
0.05 mag yields distances of 12.52 ± 0.38 and 12.11 kpc, 
both consistent with estimates from the literature using 
other methods (Table EJ. 

We checked t hat our distance modulus is consistent with 
the isochrones of lVandenBerg et al.l (|2006l ) for the metallic- 
ity of [Fe/H]zvv = —1-63 that we calculated in our analysis; 
these are overplotted on Fig. [3] and show that our results are 
broadly consistent with the theoretical isochrones of that 
metallicity. Although isochrone fitting to the CMD could 
also be used to estimate the age of NGC 1904, the quality 
of our CMD does not allow us to do this in this case, and 
we can only show that the CMD is consistent with the wide 
range of ages reported in the literature for this cluster, from 
10 to 18 Gyr (see Tabled]). 



6. The My- [Fe/H] relation 

Although recent theoretical m odels of the horizontal branch 
(HB) (e.g. IVanden Bcrg 200(3) predict t hat the relation be - 
tween My and [Fe/H] is non- linear (e.g. lCassisi et al.lll999ft . 



it has traditionally been expressed linearly in the literature 
in the form My = a [Fe/H] + (3. In this section we compare 
our values of these parameters for NGC 1904 with other 
clusters for which Fourier decomposition of RR Lyrae was 
performed. These are plotted in Fig. [HI and the data are 
listed in Table [TUJ Note that when values in the literature 
were not given on the ZW scale, we converted them to the 
ZW scale to ensure a homogeneous sample. Similarly, the 
values of My for each sample were all converted to values 
on the "long" distance scale. In several cases, this meant 
adjusting the value of My by subtracting 0. 2 to the quoted 
value mag to b ring it from the scale used bv lKovacsl (|1998|) 



to the scale of ICacciari et al.l (|2005t ). There were some ex- 
ceptions, notably for NGC 6388 and NGC 6441, as i t is un- 
clear whether the shift suggested by ICacciari et al.l (|2005l) 
applies to metal-rich clusters. The scale used in the original 
published value for each cluster is noted in Table [TO] 

Our best linear fit is My = (0.16 ± 0. 01) [Fe/H] zw + 
0.85 ± 0.02), or, using the UVES scale of ICarretta et all 
2009h . My = (0.14±0.01) [ Fe/HlTTVF,K + (0.81±0. 2). This 
is in good agreement with lArellano Ferro et al.l (|2008bl) 
who found My = (0.18 ± 0.03) [Fe/H] zw + (0.85 ± 0.05). 
It is also in good ag reement with the value found by 
iFusi Pecci et all (|1996l i. who analysed the CMDs of eight 
globular clusters in M31 to derive the absolute magnitude 
of the HB at the instability strip, and found My = (0.13 ± 
07) [ Fe/H] zw + (0.95 ± 0.09), although iFederici et all 
(l2012f) recently found My = (0.25±0.02) [Fe/H] zw +(0.89± 
0.03) for the M31 clusters, a steeper slope than our value. 
We can also compare our values to some of the many esti- 
mates of the slope a an d the intercept J3 in the literature, 
e.g. (a,0) =(0.22 0.89) (lOratton et al.ll2003H (U214 0.88 ) 
()Clementini et al.ll2003f ). Y0.18. 0.90) (jCarretta et~aT]|2000D . 
all in excellent agreement with our derived value. 

The point corresponding RR0 stars on Fig. [^appears to 
be an outlier in the distribution of My vs. [Fe/H]. Although 
this could be the result of poor light curve decomposition, 
we verified in Sec. 0] that the number of harmonics we fit 
has little influence on the parameters we derive for the RR 
Lyrae stars, and therefore on the cluster parameters, so 
this is unlikely to be the reason for the position of that 
data point on Fig. [9l Since the sample size is small, and 
error bars are quite large, we refrain from claiming that 
NGC 1904 is a clear outlier of the distribution, being only 
~ 2<T away from our best linear fit of My vs. [Fe/H]. 



7. The distribution of variables in NGC 1904 

From the finding chart of variable stars in NGC 1904 on 
Fig. we can see that the variables in this cluster are 
distributed peculiarly, with the variables lying along the 
Northwest - Southeast axis, rather than distributed ran- 
domly across the cluster as would be expected. In order to 
assess the statistical significance of this, we first compared 
the distribution of the RR Lyrae variables we detect to the 
distribution of HB stars in the cluster. This is shown in 
Fig. [lOj and the distribution of HB stars can be consid- 
ered to be spherically symmetric, unlike the distribution of 
variable stars. 

In order to assess the significance of this, we drew 10 5 
random samples of 10 stars from the HB star population in 
our data and fitted a linear function y — ax + b to each of 
these, minimising the total square perpendicular distance 
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Reference 



[Fe/Hj 2 



Method 



This work 
Kraft fc Ivans (2003) 
Carretta fe Gratton (1997) 
Kravtsov et al. (1997) 
Francois (1991) 
Gratton fc Ortolani (1989) 
Brodie fc Hanes (1986) 
Smith (1984) 
Nelles fc Seggewiss ( 1984) 
Zinn fc West (1984) 
Zinn (1980) 



-1.63 ± 0.14 Fourier decomposition of RR Lyrae light curves 

-1.68 EW of Fe II spectral lines 

-1.36 ± 0.09 Spectroscopy of red giants 

-1.76 ± 0.20 CMD analysis 

-1.46 ± 0.15 Q39 spectral index 

-1.42 ± 0.23 Absorption line strength indices 

-1.70 ± 0.23 Absorption line strength indices 
-1.43 ± 0.20 Corrected Q39 index 

-1.78 ± 0.25 Calibration of VBLUW indices 
-1.69 ± 0.09 Q39 spectral index 

-1.76 ± 0.05 Q39 spectral index 



Table 7. Different metallicity estimates for NGC 1904 in the literature. 



Reference 



fio [mag] Distance [kpc] 



Method 



This work 
This work 
Zoccali fc Piotto (2000) 
Ferraro et al. (1999) 
Kravtsov et al. (1997) 
Harris (1996) 
Zinn (1980) 
Rosino (1952) 



15.63 ± 0.06 13.35 ± 0.35 Fourier decomposition of RR0 light curves 



15.54 
15.52 
15.63 
15.57" 

15.62 Q 
16.25 



12.87 Fourier decomposition of RR1 light curves 

12.53 Theoretical luminosity function 

13.37 Magnitude of the horizontal branch 

13.00 CMD analaysis 

12.9 Globular cluster catalogue 

13.29 Q39 spectral index 

17.78 Median magnitude of RR Lyrae 



Table 8. Modulus and distance estimates for NGC 1904 in the literature. "Assuming the values E(B — V) 
Rv = 3.1. 



0.01 and 



Reference 



Kolcva et al. (2008) 
Koleva et al. (2008) 
De Angeli et al. (2005) 
De Angeli et al. (2005) 
Santos fc Piatti (2004) 
Zoccali fc Piotto (2000) 
Kravtsov et al. (1997) 
Kravtsov et al. (1997) 
Alcaino et al. ( 1994) 



Age [Gyr 



Method 



11.99 ± 0.57 Spectrum fitting 

14.02 ± 0.25 Spectrum fitting with hot stars 

10.19 ± 0.08 Age- ZW metallicity fitting 

9.92 ± 0.08 Age- CG metallicity fitting 

13.2 ± 1.1 Age-metallicity fitting 

14 RGB luminosity function 

16 ± 1 CMD isochrone fitting 

18 ± 1 CMD isochrone fitting 

16 CMD isochrone fitting 



Table 9. Age estimates for NGC 1904 in the literature. 



(TSPD) as a measure of how aligned the stars in the sample 
are. The TSPD is calculated as 

Tspp = f py^iV, (16) 

JrJ V Va 2 + 1 / 

where N = 10 is the number of variables in the sample. 
We then compared this to the same statistic for the fit 
to the RR Lyrae stars in the cluster, TSPD var . We find 
that the probability that a random draw from the overall 
distribution is equal to or lower than TSPD var is ~ 2.8%. 
The resulting distribution of TSPD is shown on Fig. QTJ 

8. Conclusions 

Using difference image analysis, we were able to obtain good 
photometry for the stars in NGC 1904, even in the crowded 
central region, where other methods struggle to overcome 
problems caused by the crowded field. Although the pho- 
tometry we have obtained for this cluster is not as accu- 
rate as that whi ch we obtained for other clusters in previ- 
ous s tudies (e.g. lArellano Ferro et al.ll2010t iBramich et al.l 
1 2 llh . we discovered a new RR1 variable and verified that 



one object is not in fact variable. Furthermore, the long 
time baseline of almost 8 years allowed us to derive pre- 
cise periods for the RR Lyrae. This in turn has enabled us 
to estimate some of the properties of the cluster through 
Fourier decomposition of the RR Lyrae light curves. 

Using this we found a metallicity for NGC 1904 of 
[Fe/HJ zw = — 1.63±0.14, or, on the scale of lCarretta et al.l 
(|2009t ). [Fe/H] UV ES = -1-57 ± 0.18. We also find distance 
moduli of /i = 15.64 ± 0.07 and jU = 15.54 for RR0 and 
RR1 variables, translating into distances of 13.26 ± 0.41 
kpc (using RR0 variables) or 12.87 kpc (using the one RR1 
variable with a good period estimate). 

Finally, we also used our CMD to check that the metal- 
licity and distance modulus we derived for NGC 1904 
is broadly consistent w ith the theoretical isochrones of 
IVandenBerg et al] (|2006l ). and found best-fit isochrones in 
agreement with the spread of ages reported in the litera- 
ture. 
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[Fc/HJzw 
(RRO) 



[Fe/Hjzw 
(RR1) 



NGC # Messier # 



M v 
(RRO) 



M v 
(RR1) 



Reference 



Distance 
scale 



NGC 1904 
NGC 6981 
NGC 5024 
NGC 5286 
NGC 5053 
NGC 6266 
NGC 5466 
NGC 6366 
NGC 7089 
NGC 7078 
NGC 5272 
NGC 4147 
NGC 6388 
NGC 6441 
NGC 6362 
NGC 6934 
NGC 5904 
NGC 6333 
NGC 6809 



M 79 
M 72 
M 53 



M 62 



M 2 
M 15 
M 3 



M 5 
M 9 
M 55 



0.41 ± 0.07 
0.62 ± 0.00 
0.46 ±0.05 
0.52 ±0.04 
0.49 ± 0.06 
0.63 ±0.03 
0.52 ±0.11 
0.66 ±0.04 
0.51 ±0.10 
0.47 ±0.03 
0.60 ± 0.02 
0.60 ±0.06 
0.66 ±0.14 
0.68 ±0.03 
0.66 ±0.01 
0.61 ±0.01 
0.61 ±0.01 
0.48 ±0.03 
0.53 ±0.09 



-1.60 ±0.15 
-1.48 ±0.03 
-1.72 ±0.06 
-1.68 ±0.15 
-1.76 ±0.13 
1.31 ±0.11* 
-1.81 ±0.12 
-0.87 ±0.14 
-1.64 ±0.11* 
-1.92 ±0.17* 
-1.59 ±0.08* 
-1.47 ±0.22* 
-1.40 ±0.16 
-1.30 ±0.13 
-1.26 ±0.03* 
-1.53 ±0.04 
-1.47 ±0.01* 
-1.85 ±0.06* 
-1.61 ±0.20* 



0.58 ±0.01 
0.57 ±0.01 
0.49 ± 0.05 
0.57 ±0.04 
0.55 ± 0.05 
0.51 ±0.03 
0.53 ± 0.06 

0.51 ±0.06 
0.52 ± 0.03 
0.57 ±0.04 
0.52 ±0.07 
0.82 ± 0.06 
0.79 ± 0.03 



-1.73 ±0.06 
-1.66 ±0.15 
-1.84 ±0.13 
-1.71 ±0.23 
-1.97 ±0.18 
-1.23 ±0.09 
-1.92 ±0.21 

-1.74 ±0.19* 
-2.10 ± 0.11* 



-0.74 ±0.23 
-0.80 ± 0.24 



This work a 

Bramich et al. (2011 ') a 

i AreHano i l^ejTo 1 ^LiiL (2011) a 

Zorotovic et al. (2010) b 

^rehano i l^erro i ^£ - aL i (2010) a 

Contreras et al. (2010) b 

Arellano Ferro et al. (2008b) a 

Arellano Ferro et al. (2008a) c 

Jjazaro_et ii aL i (2006) b 

Arellano Ferro et al. (2006) b 

Cacciari et al. (2005 ) a, d 

Arellano Ferro et al. (2004) b 

Pritzl et al. (2002) 6 1 

Pritzl et al. (2001) 6 1 

~ (2001) b 
b 



^lecl^e^a^ 

Kaluznv et al. (2001) 
Kaluznv et al. (2000) b 
Clement fc Shelton (1999a) b 
Olech et al. (1999) b 



Table 10. My and [Fe/H]zw estimates for other clusters in the literature for which these parameters were calculated 
through Fourier decompo sition of the light curv es of RR Lyrae stars. All values of My were converted to a value consistent 
with t he distance sca l e of | Cacciari et al.l (l2005|) i f another scale was us ed. Distance scales used in the original publication 
are: a lCacciari et al.l (120051 ): frlKovacsl (119981 ) : c | Cacciari et al.l (|2005|) but with a shift of 0.18 mag rather than 0.2 mag 
compa red to lKovaca (jl998l ): rf lCacciari et al.l (|2005l ) but with a shift of 0.23 mag rather than 0.2 mag compared to lKovaca 
(1 19981). 1 The values of M y for these two clusters were not shifted because of uncertainty as to whether the shi ft given 
by Cacciari et all (|2005l ) applies to metal-rich clusters. *those metallicity values were published on the scale of iJurcsikl 
(119951) : values listed here are on the ZW scale and were converted using Eq. (JT]) . 




1 .5 
[Fe/H] 



Figure 9. < My > vs [Fc/H] for clusters in the litera- 
ture for which Fourier decomposition of the RR Lyrae vari- 
ables was performed. Points obtained by analysing the light 
curves of RRO variables are shown as black filled circles, 
while points from analysis of RR1 variables are marked as 
red filled circles; the data for NGC 1904 are shown as a 
green cross (RRO variables) and a green filled circle (RR1 
variables). The best linear fit is shown as a blue line, with 
coefficients given in the text. All points are plotted with 
l-cr error bars. 



-24:32:15 - 



*:27:20 L 

05:23:38 



Figure 10. Plot of the spatial distribution of all HB stars 
(black filled circles) and RR Lyrae stars (red filled circles) 
on our images. North is up and East to the right. The direc- 
tions to the Galactic Centre (GC) and Canis Major dwarf 

;alaxy (CMa) are al so shown. The tidal radius r t ~ 500" 

Lanzoni et al.l 120071 ) lies outside our plot. 
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A x 1 CT 6 '-■ 




2x10 5 4x1G 5 6x'0 5 8x10 5 1 x10 6 

Tc;a sqjare perpendicular distance (TSPD) 



Figure 11. Histogram distribution of the total square per- 
pendicular distance (TSPD) for a linear fit to 10 5 random 
samples of 10 HB stars, compared to TSPD var , the TSPD 
of the fit to the RR Lyrae stars in NGC 1904 (red vertical 
dotted line) . The cumulative probability distribution is also 
plotted as a dashed line; the probability that the value of 
TSPD is equal to or lower than TSPD var is - 2.8%. 
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